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In this paper, we show how to use low-fidelity operations to control the dynamics of quantum 
systems. Noisy operations usually drive a system to evolve into a mixed state and damage the 
coherence. Sometimes frequent noisy operations result in the coherent evolution of a subsystem, 
and the dynamics of the subsystem is controlled by tuning noisy operations. Based on this, we find 
that universal quantum computation can be carried out by low-fidelity (fidelity < 90%) operations. 
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I. INTRODUCTION 

Manipulating quantum systems coherently is impor¬ 
tant for quantum computation [1], which is believed to 
have non-trivial advantages over classical computation. 
In recent years, persistent quantum memories (e.g. [2]) 
and precise quantum operations have been demonstrated 
with individual qubits or clusters of a few qubits [3-5]. 
However, quantum computation is still a challenge, and 
one of the main obstacles is the difficulty to maintain 
the fidelity of quantum operations when many qubits are 
assembled together. Some alternative models other than 
the standard model of quantum computation have been 
proposed to exploit different mechanisms of processing 
quantum information, e.g. adiabatic quantum computa¬ 
tion [6], measurement-based quantum computation [7], 
and dissipation quantum computation [8]. In this paper, 
we will show a protocol of quantum computation utiliz¬ 
ing noisy operations, i.e. the fidelity of operations is lower 
than 90%. 

Decoherence and operation imperfections always in¬ 
duce some errors on qubits. The number of errors in¬ 
creases with the time, the number of operations, and the 
number of qubits without error correction, which could 
finally cause failures of quantum computing. The theory 
of fault-tolerant quantum computation (FTQC) predicts 
a threshold of the error rate: if errors occur with a rate 
below the threshold, errors are correctable and the com¬ 
putation is reliable [9]. For topological codes, the error 
rate threshold is about 1% (one error in a hundred op¬ 
erations) [10, 11], which is among the best records of 
the threshold. By combining the idea of noisy-operation 
quantum computation proposed in this paper and a topo¬ 
logical code, we find that a high error rate > 10% is tol¬ 
erable for realistic frequency of operations and coherence 
time. 

In our proposal of quantum computing with noisy op¬ 
erations, we consider a system composed of two coupled 
subsystems A and Q, in which the subsystem A (actua¬ 
tor) directly suffers noisy operations, and the subsystem 
Q is a register storing the quantum state for processing. 
By frequently performing noisy operations, the subsys¬ 
tem A is decoupled from the subsystem Q and fixed in 
a mixed state. In this case, the subsystem Q evolves 


solely and coherently, however, its dynamics depends on 
the fixed state of the subsystem A [12]. Therefore, by 
changing noisy operations to alter the fixed state of the 
subsystem A, one can effectively tune the dynamics of 
the subsystem Q. Based on this new idea of the control 
of a quantum system, we will show that universal and 
scalable quantum computing can be achieved with noisy 
operations. By investing the fault-tolerance thresholds, 
we hnd that these operations can be very noisy. 

The idea of indirect control using projective measure¬ 
ments [13], completely-controlled dynamics [14, 15], or 
initializations [16] of an ancillary system has been stud¬ 
ied theoretically and has applications in hybrid systems 
composed of electron spins and nuclear spins [17-21]. 
If those frequent operations performed on the actuator 
subsystem are unitary rather than noisy, both subsys¬ 
tems can evolve independently and coherently, which is 
known as the dynamical decoupling [22]. In this paper, 
we first propose to use noisy operations with low fideli¬ 
ties rather than high-quality quantum operations as the 
resource for processing quantum information. In this sce¬ 
nario, the actuator subsystem is always fixed in a mixed 
state independent of the register subsystem. 


II. QUANTUM GATES BASED ON NOISY 
OPERATIONS 

We suppose that the free time evolution of the system 
is given by the Hamiltonian H, in which two subsystems 
are coupled with each other. The set of operations that 
can be performed on the actuator subsystem includes the 
initialization operation !(•) = X[j|0)(jU ’ |j)(0U 
unitary operations {U{-) = U ■ C/'^}. Here, {Ij)^} are 
the basis states, and {U} are unitary operators of the 
actuator subsystem. With imperfections, operations ac¬ 
tually performed, which are respectively denoted by T' 
and {W} for the initialization and unitary operations, 
are different from those ideal operations I and {U}. An 
imperfect initialization prepares the actuator in a mixed 
state rather than the pure state |0)^, and an imperfect 
unitary operation corrupts as well as rotates the state of 
the actuator subsystem. 

When operations are ideal, by frequently repeating the 
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initialization X followed by a unitary operation the ac¬ 
tuator subsystem is frozen in the state \'4’)A = U\Q) a- Be¬ 
cause only the actuator subsystem is completely frozen, 
the whole system evolves in a subspace with an effective 
Hamiltonian predicted by 

the quantum Zeno effect theory [23, 24]. Here, the com¬ 
bination of the initialization and the unitary operation 
is similar to a projective measurement with \'tp )a as the 
output state [29]. When operations are imperfect, by fre¬ 
quently repeating the noisy initialization I' followed by 
a noisy unitary operation lA', the actuator subsystem is 
frozen in a mixed state pu = lA'X' {^a/ dA) rather than 
a pure state, where 1a is the identity of the actuator 
subsystem, and Aa is the dimension of the Hilbert space. 
Here, the combined operation lA'I' is a projector in the 
operator space {U'X'U'T = lA'I'), which leads to the 
operator quantum Zeno effect [12]. Then the whole sys¬ 
tem evolves, similar to the normal quantum Zeno effect, 
with an effective Hamiltonian Hyi{U) = 1a®Xt:a{puH)- 
In this effective Hamiltonian, two subsystems are decou¬ 
pled, and the effective Hamiltonian of the register sub¬ 
system is given by Hq{U) = Tta{puH)- Therefore, by 
changing the noisy unitary operation IA\ i.e. changing the 
fixed state pu of the actuator subsystem, we can control 
the dynamics of the register subsystem. We would like 
to remark that, in order to precisely control the regis¬ 
ter subsystem, both the system Hamiltonian and noise 
in operations, i.e. the map between the operation U and 
the effective Hamiltonian Hq{U), must be known. 

In the standard model of quantum computation, a uni¬ 
versal set of quantum gates includes a set of single-qubit 
gates and at least one two-qubit entangling gate such as 
the controlled-phase gate [1]. In the following, we will 
give some examples of implementing gate operations on 
the register subsystem with noisy operations on the ac¬ 
tuator subsystem. In Sec. HIB, we will show how to 
prepare and measure the register subsystem. 


A. Single-qubit gate 

As the first example, we consider two qubits coupled 
via the Heisenberg interaction = J{a\aQ -I- -I- 

^a^q)- These two qubits are the actuator subsystem 
A and the register subsystem Q, respectively. We model 
imperfect operations as a combination of ideal operations 
and depolarizing noise. An operation with depolarizing 
noise reads O' = where O is the ideal operation, 
and £e(-) = (1 - 3e/4) • -k(e/4)(cr"^ • -k cr« • -f cr^ • a^) 
is the depolarizing operation. In this model of noise, 
the ideal operation is performed with the probability 
1 — e, while the state of the qubit is completely de¬ 
stroyed, i.e. turns into the maximally mixed state, with 
the probability e. Under frequent noise operations, the 
fixed state of the actuator qubit can always be writ¬ 
ten as Pu = (1/2)1^ + PxC^A + + Vz'^a- Then, 

the corresponding effective Hamiltonian of the register 
qubit is Hq(U) = JijpxCrq + Py(JQ -k p^crg). If only 



FIG. 1: Error rate of a single-qubit phase gate cr^ on the 
register qubit implemented with frequent noisy operations on 
the actuator qubit. F is the entanglement fidelity [25] of the 
phase gate. The error rate (1 — F) always increases with the 
depolarizing rate and decreases with the frequency of noisy 
operations. Three contour lines correspond to phase-gate fi¬ 
delities 99%, 99.9%, and 99.99%, respectively. See Sec. A for 
details. 


the initialization is performed, the fixed state is p^ = 
(1 — ei/2)j0)(0lA + ei/2ll)(llA, and the corresponding ef¬ 
fective Hamiltonian is Hq{1) = (1 — ei)JaQ, where is 
the depolarizing rate of the noisy initialization. By com¬ 
bining the initialization with a Hadamard gate, the fixed 
state is changed to pn = {2- ei- eh + eie;i)/2l-k)(-kU + 
(ci -k — eieh)/2\—){—\A, and the corresponding effec¬ 
tive Hamiltonian is Hq(H) = (1 — ei)(l — eh)JcrQ. Here, 

1±)^ = (I/-\/2)(j0)yi ± 1 I)a), and eh is the depolarizing 
rate of the noisy Hadamard gate. Therefore, in this ex¬ 
ample, the initialization and Hadamard gate on the ac¬ 
tuator qubit are enough for universal single-qubit gates 
on the register qubit. 

When the frequency of noisy operations is finite, the 
evolution of the register subsystem is not exactly de¬ 
scribed by the effective Hamiltonian, which may result in 
errors in the register subsystem. In the previous example 
of implementing single-qubit gates, a phase gate CTq on 
the register qubit is equivalent to the time evolution de¬ 
scribed by exTp[—iHQ{l)t/h] with t = Ii7r/[2(1 — e^jJ]. 
When the frequency of noisy operations is finite, the 
fidelity of the phase gate decreases with the depolariz¬ 
ing rate and increases with the frequency of noisy op¬ 
erations as shown in Fig. 1. We find that even if the 
depolarizing rate is very high (e.g. 80%), a high-fidelity 
(e.g. 99%) phase gate can still be achieved (e.g. with the 
frequency ~ 1.5 x W^J/h). Actually, when the frequency 
approaches infinite, the depolarizing rate only affects the 
time cost of implementing the phase gate, i.e. is limited 
by the coherence time. 
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FIG. 2: Networks for universal quantum computation with 
noisy operations, (a) A network that allows universal quan¬ 
tum computation. Each register qubit (round) is coupled with 
one actuator qubit (triangle) for controlling single-qubit gates 
and four actuator qubits (squares) for controlling two-qubit 
gates, (b) A network with only two-body interactions. Each 
double-round register represents a four-state system encoding 
two qubits. The triangle actuator is coupled with only one 
qubit of the register (the core qubit), and the square actua¬ 
tor is coupled with both qubits (the core qubit and the port 
qubit). Double-round registers are connected by two-body 
interactions with actuator qubits (rings) for controlling the 
state transfer. 


B. Two-qubit gate 

For the two-qubit entangling gate, we consider three 
qubits coupled via a three-qubit Ising interaction Hi = 
J^A^Qi^Q 2 ^ where qubits Q1 and <52 form the register 
subsystem, and qubit A forms the actuator subsystem. 
With frequently initializing the actuator qubit, the dy¬ 
namics of two register qubits is given by the effective 
two-qubit Ising interaction Hqq(1) = J(1 — ei)o'QiCrQ 2 - 
And the time evolution expl—iHQQ(l)t/h] with t = 
/i7r/[4(l — ei)J] gives the two-qubit phase gate Rzz = 
(1 — ZCTQ]^tTQ 2 )/-\/ 2 ) which can maximally entangle two 
register qubits and is identical to the standard controlled- 
phase gate up to single-qubit phase gates. In Sec. Ill A, 
we will show how to implement a two-qubit gate with 
only two-body interactions. 


III. UNIVERSAL QUANTUM COMPUTATION 
WITH NOISY OPERATIONS 


To have a universal quantum computer, we need to in¬ 
tegrate single-qubit gates and two-qubit gates in the same 
scalable network. For example, in the two-dimensional 
network shown in Fig. 2 (a), each register qubit (round) is 
coupled with an actuator qubit (triangle) via the Heisen¬ 
berg interaction Hu for single-qubit gates, and each pair 
of neighbouring register qubits are coupled with an ad¬ 
ditional actuator qubit (square) via the three-qubit Ising 
interaction Hi for two-qubit gates. In this network, all 
register qubits form the register subsystem, and all ac¬ 
tuator qubits form the actuator subsystem. We have 
discussed how to implement single-qubit gates and the 
two-qubit phase gate via the Heisenberg interaction and 
the three-qubit Ising interaction, respectively. Gate op¬ 


erations on register qubits are switched off (switched to 
1) by frequently performing the single-qubit twirling op¬ 
eration £i (a depolarizing operation with the rate 1) on 
corresponding actuator qubits. The twirling operation is 
equivalent to randomly performing Pauli gates. With fre¬ 
quent twirling operations, an actuator qubit is frozen in 
the maximally-mixed state p = (|0)(0|^-|-|1)(1|a)/2, and 
then the corresponding effective dynamics is switched off 
(the effective Hamiltonian oc 1). By changing between 
twirling operations and other noisy operations, i.e. ini¬ 
tialization and unitary-gate operations, a quantum cir¬ 
cuit can be implemented in this network. 


A. Two-body-interaction model 


The network for universal quantum computation can 
also be built with only two-body interactions. In the 
network shown in Fig. 2 (b), the elementary unit is a 
complex composed of two actuator qubits (triangle and 
square) and a four-state register particle (double round), 
e.g. a spin-3/2 particle. Here two qubits are encoded 
in each register particle, which are respectively called 
the core qubit and the port qubit. The triangle actu¬ 
ator qubit is coupled to the core qubit via H^. And 
the square actuator qubit is coupled to both qubits via 
Hi, which, however, is a two-body interaction because 
two qubits are in the same particle. Register parti¬ 
cles are connected via the third kind of actuator qubits 
(rings). The coupling between a pair of neighbouring 
register particles and a ring actuator qubit is the XY- 
interaction H^y = + ^P 2 ) + + ^P 2 )]j 

where complex-1 and complex-2 are two neighbouring 
complexes. Cl and P2 are respectively the core qubit 
of complex-1 and the port qubit of complex-2, and A is 
the ring actuator qubit. In this network, the informa¬ 
tion processed in the quantum computing is stored in 
core qubits. Single-qubit gates on core qubits are im¬ 
plemented via triangle actuator qubits. Two-qubit gates 
are achieved by the state transfer between core qubits 
and port qubits in neighbouring complexes. The XY- 
interaction Hxy drives a time evolution ending up with 
a swap gate between two qubits Cl and P2 at the time 
t = Htt/ 2^f/2J. By frequently performing the twirling op¬ 
eration £i on the ring actuator qubit, the XY-interaction 
can be effectively switched off, and complexes are decou¬ 
pled. When a two-qubit phase gate on core qubits Cl 
and C2 is required, the twirling operation is turned off 
for the time t. Then the free time evolution transfers 
the state from the core qubit Cl to the port qubit P2. 
After a local two-qubit phase gate on qubits C2 and P2 
(via the square actuator qubit) and another swap gate, a 
two-qubit phase gate on Cl and C2 is achieved [see the 
circuit in Figs. 3(a) and 3(b)]. We would like to remark 
that, each gate given by the time evolution 

e _ -(cTaCTci + ^A^P2 + ^Cl^P2 ~ 1)SWAP 
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FIG. 3: Circuits for (a) the effective two-qubit phase gate 
on neighbouring complexes, (b) the corresponding operations 
on actuator qubits, and (c) the measurement distillation. In 
(b), zigzag lines represent frequent twirling operations 1 
denotes turning off the twirling operation, and T' denotes fre¬ 
quent initialization operations. In (c), both the data qubit 
and the ancillary qubit are register qubits, and controlled 
phase gates are implemented with frequent noisy operations 
on the actuator subsystem. In each round of the distillation, 
the ancillary qubit is initialized in the state |-|-) and measured 
in the basis. 


has an additional phase to the net swap gate SWAP = 
i^Ci^P 2 + ®’ci^P 2 + ^ci^P 2 + l)/2 depending on initial 
states of the port qubit P2 and the ring actuator qubit 
A. Fortunately, the additional phases attached with two 
swap gates are cancelled with each other as they commute 
with the two-qubit phase gate. Therefore, the overall 
operation on two core qubits Cl and C2 is independent 
of the initial states of the port qubit P2 and the ring 
actuator qubit A. 


B. Initialisation and measurement of register 
qubits 


Besides gate operations, preparation and readout of 
the state of qubits are also required by quantum com¬ 
puting. Here we suppose that register qubits can be di¬ 
rectly initialized and measured but with a low fidelity. 
If gate operations on register qubits can be implemented 
with high fidelity, we can measure the state of a regis¬ 
ter qubit precisely with the help of a distillation circuit 
[see Fig. 3 (c)]. In the ideal situation that there is no 
error in gate operations on register qubits, an accurate 
effective measurement on the data qubit can always be 
achieved by repeating the circuit for many times, and the 
state of the data qubit can be read (in the basis) from 
the majority of measurement outcomes of the ancillary 
qubit. When gate operations also have errors, the proto¬ 
col of reading the data qubit can be adapted to improve 
its fidelity (see Sec. A). Once the state of the data qubit 
is successfully measured, the data qubit is initialized in 
the state either |0) or |1). 



FIG. 4: Thresholds of the depolarizing rate in noisy- 

operation quantum computation. From bottom to top the 
lines correspond to distilling measurements for 9, 11, 13, 15, 
and 17 rounds, respectively. Below the threshold, errors on 
register qubits can be corrected with the error correction code. 
We would like to remark that thresholds could be further im¬ 
proved by considering higher frequencies and more rounds of 
distillations. To obtain these thresholds, we have assumed 
that operations performed on actuator qubits and register 
qubits are all noisy, and noises are depolarized with the same 
depolarizing rate. 


IV. FAULT-TOLERANT QUANTUM 
COMPUTATION AND THRESHOLDS 


Finally, we would like to show FTQC thresholds of 
noisy operations in Fig. 4. For these thresholds, we have 
considered using a three-dimensional version of the net¬ 
work shown in Fig. 2 (a) to generate a topology-protected 
cluster state (see Sec. A for details), in which errors can 
be corrected if less than 3% of qubits are affected by er¬ 
rors [26]. We also have assumed that all operations are 
noisy, and noises are depolarized with the same depo¬ 
larizing rate. In the result, for a frequency of 10^J//i, 
the depolarizing rate threshold is about 20.1%, which 
corresponds to the error rate 10.05% for the noisy ini¬ 
tialization and measurement and the error rate 15.075% 
for noisy unitary operations. To obtain the thresholds, 
we have neglected the environment-induced decoherence. 
In the 17-round distillation case, all operations on regis¬ 
ter qubits from initialization to measurement are finished 
within the time 7{J/h)~^. Hence, if the coherence time is 
much longer than 7{Jjh)~^j3% ~ 233(J//i)“^, errors in¬ 
duced by decoherence occur with a rate much lower than 
3% and only reduce the thresholds slightly. We have also 
neglected the fluctuation of interactions and operation 
noises. A small fluctuation (<C 1%) of the interaction 
strength J and the noise parameter e will not affect the 
threshold significantly. 
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V. DISCUSSIONS AND SUMMARY 

A candidate for realizing noisy-operation quantum 
computation is the kind of hybrid systems composed of 
electron spins and nuclear spins [17-21]. Electron spins 
in quantum dots can be initialized and manipulated in 
tens of picoseconds [27, 28] and can play the role of ac¬ 
tuators. Nuclear spins are usually well decoupled from 
the microwave and optical pumping for electron-spin op¬ 
erations and can have an extremely long coherence time 
(e.g. hours [2]), hence, can play the role of registers. With 
the coupling strength of 1 MHz between electron spins 
and nuclear spins, provided noisy operations on electron 
spins can be performed in 100 ps, one can operate nuclear 
spins in ^ 1 /rs. 

In summary, we have discussed how to use noisy op¬ 
erations to process quantum information and obtained 
FTQC thresholds in an example model. We find that 
fidelities of operations can be even lower than 90%. Our 
results provide a way to achieve quantum computation 
by boosting the operation frequency rather than the op¬ 
eration fidelity. However, we would like to remark that 
this protocol of noisy-operation quantum computation 
cannot replace error correction codes. Actually, the er¬ 
ror correction is an important component of the over¬ 
all protocol. Because quantum computation is the most 
complicated task among other applications of quantum 
technologies, we believe that the same idea can also be 
used in quantum communication and quantum sensing. 
In this paper, we have focused on isotropic Heisenberg in¬ 
teractions and Ising-type interactions as examples. The 
same idea can be applied to other types of interactions, 
e.g. interactions between electron spins and nuclear spins 
could be anisotropic. Although we have only discussed 
depolarising noise in detail, noises are not restricted to a 
specific type. Necessary conditions for noises are given in 
Sec. B, and only a small set of noises, in which operations 
are extremely noisy, are not suitable for noisy-operation 
quantum computation. 


FIG. 5: A three-dimensional NOQC network for building the 
topology-protected cluster state, (a) An elementary cubic of 
the topology-protected cluster state. On the cluster state, 
each vertex represents a qubit initialized in the state |-|-), and 
each edge represents a controlled-phase gate on corresponding 
qubits, (b) Each qubit on the cluster state corresponds to a 
data qubit (d.q.) on the NOQC network. An ancillary qubit 
(a.q.) is associated with every data qubit for initialization and 
measurement distillations, (c) Neighbouring data qubits are 
connected according to the cluster state lattice. Here, both 
data qubits and ancillary qubits are register qubits, which are 
controlled by noisy operations applied on actuator qubits. 


To build the topology-protected cluster state, we need 
a three-dimensional version of the network shown in 
Fig. 2 (a). On the three-dimensional network (see Fig. 5), 
each cluster-state qubit corresponds to a data register 
qubit, and an ancillary register qubit is attached to ev¬ 
ery data qubit for initialization and measurement distil¬ 
lations. Each register qubit is coupled with a triangle 
actuator qubit for single-qubit gates, and each pair of 
connected register qubits are coupled with a square actu¬ 
ator qubit for two-qubit gates. The interactions between 
qubits, iJn and Hj, are discussed in Sec. H. We suppose 
that, the initialization and measurement in the tr^ basis 
{[0), Jl)}, the Hadamard gate, and Pauli gates could be 
applied on actuator qubits, the initialization and mea¬ 
surement in the tr^ basis {[+), |—)} could be applied on 
register qubits, and all of these operations have the same 
depolarizing rate e and could be repeated with the same 
frequency /. Here, a noisy measurement is modelled as 
a depolarizing error followed by a perfect measurement. 
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Appendix A: Threshold simulation 
1. Model 

To obtain the thresholds of fault-tolerant quantum 
computation (FTQC) on the noisy-operation quantum 
computation (NOQC) architecture, we consider build¬ 
ing the topology-protected cluster state [see Fig. 5 (a)], 
which can tolerate phase errors with a rate of 3% [26]. 


2. Circuit 

The cluster state is prepared and measured with the 
overall circuit shown in Fig. 6. First, data qubits are ini¬ 
tialized in the cr^ basis by repeating the distillation cir¬ 
cuit. After initialization, data qubits are in states either 
jo) or jl) depending on measurement outcomes in distil¬ 
lation circuits. Flip operations could be used to align all 
data qubits to the state JO), but they are not necessary 
because one can update the basis rather than physically 
implement flip operations. Second, data qubits are fur¬ 
ther prepared in the state j-f) via single-qubit gates, and 
then a cluster state can be built on the network with two- 
qubit phase gates. Finally, data qubits are measured in 
the basis by another set of distillation circuits. In 
Fig. 6, qubits on the prepared cluster state are effec- 
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tively measured in the cr^ basis because of single-qubit 
gates between two-qubit phase gates and measurement 
distillation circuits. 

For each round of the distillation (each blue box in 
Fig. 6), ancillary qubits are initialized in the state |-|-) 
and measured in the basis. The two-qubit gate 
R'zz = Rz.., Rzz = (1 - aiq.)/2 - *(a + 
is a phase gate on the ancillary qubit depending on the 
state of the data qubit and equivalent to the standard 
controlled-phase gate upto a single-qubit phase gate on 
the data qubit. Here, cr| ^ and cr^.q. Pauli opera¬ 
tors of the data qubit and the ancillary qubit, respec¬ 
tively, Rz = (1 — is a single-qubit phase gate, 

which is implemented by frequently performing noisy ini¬ 
tializations on the corresponding triangle actuator, and 
the two-qubit phase gate Rzz is implemented by fre¬ 
quently performing noisy initializations on the corre¬ 
sponding square actuator. As Rz^,^ and Rzz commute 
with each other, the two gates are implemented at the 
same time. Each distillation circuit is effectively a mea¬ 
surement of the data qubit in the basis, hence by 
repeating the circuit, data qubits are initialized or mea¬ 
sured in the basis. 

The single-qubit gate RzRx could be used to ro¬ 
tate the data-qubit state from |0) (|1)) to |-|-) (|—)) for 
the cluster-state generation. And the single-qubit gate 
RxR^z could be used to rotate the data-qubit state from 
|-|-) (|—)) to |1) (|0)) for effective measurements in the 
basis. Here, the gate Rx = (1 — is imple¬ 

mented by frequently and alternatively performing noisy 
initializations and Hadamard gates on the correspond¬ 
ing triangle actuator. In the circuit shown in Fig. 6, 
Rz and are cancelled by each other as they com¬ 
mute with the two-qubit phase gate Rzz- The standard 
controlled-phase gate Az = (1 —crf)/2-|-(l-|-crf)cr|/2 for 
the cluster-state generation is replaced by the two-qubit 
phase gate Rzz = (1 — icrfCTf )/-\/2. There are a total of 
four controlled-phase gates applied on each data qubit, 
which are implemented at the same time as they com¬ 
mute with each other. Because Az = R^^iRz^R^z, the 

four single-qubit phase gates are cancelled as = Ij 
and each controlled-phase gate Az is replaced by a two- 
qubit phase gate Rzz- 


3. Errors 

For preparing the cluster state and measuring cluster- 
state qubits in the basis, the noisy-operation- 
controlled operations (NOCOs) implemented on register 
qubits includes R'zz^ Rx, Rzz, and decoupling opera¬ 
tions. Because we are interested in the case that the 
rate of errors on each cluster-state qubit is ~ 3%, the 
error rate of each individual NOCO is < 3%. There¬ 
fore, the probability of two errors occurring on the same 
cluster-state qubit but induced by different operations is 
< 0.09% ^ 3%. In the numerical simulations of error 


rates, we will neglect the possibility of the case that two 
errors occur on the same cluster-state qubit, i.e. when we 
study the errors induced by one NOCO, we suppose all 
other NOCOs are perfectly performed. 

For a NOCO R implemented with the noisy operation 
U'R and the interaction H, the operation actually per¬ 
formed on register qubits reads ??.'(•) = {-(S) 

Pa)], where T(-) = denotes the free 

time evolution, and / is the frequency of repeating noisy 
operations. Here, we have supposed that the time of 
performing a noisy operation is much shorter than the 
time interval between two noisy operations. pA is the 
initial state of the actuator, which could be chosen as 
lyi/2. When the unitary operation corresponding to W 
is 17 = 1, only X' will be performed on actuator qubits. 
If Rq is the effective Hamiltonian, and the desired oper¬ 
ation R = the number of noisy operations N is 

the largest integer that does not exceed t/. It is similar 
for decoupling operations, where WX' is replaced by the 
twirling operation E\. 

In the ideal case, i.e. the frequency of noisy operations 
is infinitely high, the operation TV — 7?., where the ideal 
operation ??.(•) = When the frequency is finite, the 

actually performed operation is different from the ideal 
operation and can always be expressed as TV = ErTZ. 
Here, the superoperator Er denotes the errors in the 
NOCO i?, which is found using the Choi-Jamiolkowski 
isomorphism in our numerical simulations. We would 
like to remark that the error rates of a single-qubit phase 
gate shown in Fig. I are obtained using the method de¬ 
scribed in this section. 


4. Distillation and error correction 

The high-fidelity initialization and measurement of 
data qubits could be achieved by repeating the distil¬ 
lation circuit (blue box in Fig. 6). The initialized state 
and measurement outcome of the data qubit are read 
from measurement outcomes of the ancillary qubit. If the 
distillation circuit is repeated for n times in the measure¬ 
ment (initialization) circuit, there are n measurements of 
the ancillary qubit and a total of 2” different sets of out¬ 
comes. Corresponding to each set of n ancillary-qubit 
measurement outcomes, o, the input (output) state of 
the data qubit is |0) with the probability and |1) with 
the probability I — qg- Then, when the outcomes d occur, 
the input (output) state is likely to be |0) ii qg > 1 — qg, 
or |1) if gg < 1 — qg. Therefore, the distillation fails with 
the probability pf = '^gPg^R^{qg, 1 — qg}, where pg is 
the probability of outcomes o. 

In a FTQC algorithm of the topology-protected clus¬ 
ter state [26], most of cluster-state qubits (vacuum) are 
measured in the cr® basis, i.e. proceeded following the 
circuit shown in Fig. 6. Except for vacuum qubits, some 
cluster-state qubits (defects) are measured in the basis 
for defining the computation algorithm, and some other 
cluster-state qubits (singular qubits) are measured in the 
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FIG. 6: Overall circuit for the preparation and measurement of a cluster state on the NOQC network. In this circuit, the two 
cluster-state qubits are effectively measured in the basis, which is the most likely case in the measurement-based quantum 
computation on the topology-protected cluster state. 


basis of (cr^zhcr^)/v^ for inputting magic states. The cir¬ 
cuits for defect qubits and singular qubits are slight dif¬ 
ferent from Fig. 6, i.e. one more single-qubit gate on the 
data qubit needs to be added after the second Rx gate. 
Although errors on defect qubits and singular qubits are 
different from errors on vacuum qubits due to the addi¬ 
tional gate, the different is small and does not affect the 
threshold [26]. Therefore, the threshold is mainly deter¬ 
mined by errors generated in the circuit in Fig. 6. 

By numerically simulating the errors in the circuit in 
Fig. 6 and comparing the effective phase error rate, con¬ 
tributed by both phase errors on the cluster state and 
errors of effective data-qubit measurements in the ba¬ 
sis, with the threshold 3%, we find the FTQC thresholds 
of NOQC shown in Fig. 4. We would like to remark that, 
in our model, phase errors on the same sub-lattice [26] 
are almost independent, i.e. the correlated errors occur 
with a rate < 0.03% <C 3% near the threshold, hence the 
correlations only affect the threshold slightly and could 
be neglected. 

Appendix B: Necessary conditions for noises 

In the main text, we have only discussed depolariz¬ 
ing noise in detail. Here we will show that NOQC is 
not restricted to depolarizing noise. We take the NOQC 
network shown in Fig. 2 (a) as an example. And we 
suppose that the noise in initialization operations on tri¬ 


angle actuator qubits is £i, the noise in Hadamard gates 
on triangle actuator qubits is ^h, and the noise in ini¬ 
tialization operations on square actuator qubits is £s- 
Then, with frequent initialization operations, a triangle 
actuator qubit is frozen in the state pi = £^i(l0)(0l); with 
frequent combined (initialization -b Hadamard gate) op¬ 
erations, a triangle actuator qubit is frozen in the state 
Ph = £n{HpiH)', and with frequent initialization op¬ 
erations, a square actuator qubit is frozen in the state 
PS = 'S’s(|0)(0j). To implement universal single-qubit 
gates on register qubits, states pi and pn must be po¬ 
larized in different directions, i.e. 

Tr(crpi) X Tr(crpH) ^ 0, (Bl) 

where 

cr ^ axi + ayj + cr^k. (B2) 

To implement the two-qubit phase gate, we need 

Tr(a,ps) ^ 0. (B3) 

Equations (Bl) and (B3) are necessary conditions for the 
noise in operations. These conditions are not satisfied 
only in extreme cases, i.e. i) a Hadamard gate cannot ever 
change the polarization of a qubit initialized in the com¬ 
putational basis, or ii) after an initialization operation, 
states jo) and jl) still occur with the same probability. 
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